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Abstract 

We review results from searches for mixing and CP violation in the D Q -D° 
system. No evidence for mixing or CP violation is found, and limits are set for 
the mixing parameters x, y, x' , y' , and several CP-violating parameters. 



1 Introduction 

Despite numerous searches, mixing between D° and D° flavor eigenstates has 
not yet been observed. Within the Standard Model (SM), the short-distance 
"box" diagram (which plays a large role in K°-K° and B°-B° mixing) is 
doubly-Cabibbo-suppressed (DCS) and GIM-suppressed; since the D° decay 
width is dominated by Cabibbo-favored (CF) amplitudes, D°-D° mixing is 
expected to be a rare phenomenon. Observing mixing at a rate significantly 
above the SM expectation could indicate new physics. 



The formalism describing D°-D° mixing is given in several papers. ^' ^) 
The parameters used to characterize mixing are x — Am/T and y = Ar/(2r), 
where Am and Ar are the mass and decay width differences between the two 
mass eigenstates, and T is the mean decay width. Within the SM, x and y 
are difficult to calculate as there are long-distance contributions. For m q » 
Aq CD , these contributions can be estimated using the heavy-quark expansion; 
however, m c may not be large enough for this calculation to be reliable. Current 
theoretical predictions span a wide range: \x\ ~ \y\ ~ (10~ 7 to I0~ 2 ), with 
the majority being < ICT 3 . 

For decay times t< 1 /Am, 1 / Ar, which is well-satisfied for charm decay, 
the time-dependent D°(t)—*f and D°(t)—*f decay rates are 
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where A = (q/p)(Af /Af), A = (p/q)(Aj/Af), q and p are complex coefficients 
relating flavor eigenstates to mass eigenstates, and Af (Af) and Af(Af) are 
amplitudes for a pure D° (D°) state to decay to / and /, respectively. 

In this paper we discuss five methods used to search for D°-D° mixing 
and CP violation (CPV). These methods use the following decay modes: 1 
semileptonic D° — > K + 1~ v decays, decays to CP-eigenstates K + K~ and 7r + 7r~, 
DCS D ^K+tt- decays, D° -» K% tt+tt" decays, and multi-body DCS D° -» 
K + n(ir) decays. A newer method based on quantum correlations ^) in e+e~ — > 
V>"(3770) — > D°D° production is not discussed here. The flavor of a Z?° when 
produced is determined by requiring that it originate from a D* + — ► D°tt^ 
decay; the charge of the low momentum ("slow") 7r+ determines the charm 
flavor at t = 0. As the kinetic energy released in £)* + — > -D°7r+ decays is only 
5.8 MeV (very near threshold), requiring that Q = M K7rv — M Kn — m, be 
small greatly reduces backgrounds. 



1 Charge-conjugate modes are implicitly included throughout this paper un- 
less noted otherwise. 
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2 D°(t)^K^+£-v Semileptonic Decays 

Because the + l~v final state can only be reached from a D° decay, observ- 
ing D°(t)— >K(*} + (rv would provide clear evidence for mixing. In Eq. (1) only 
the third term is nonzero; integrating this term over all times and assuming 
\q/p\ — 1 (i.e., neglecting CPV in mixing) gives 

/ R{D° -> K+lv) dt x 2 + y 2 

j R(D° ^ K-lv) dt ~ 2 = T °' ( } 

Several experiments ^' ^) have used this method to constrain ; the 
most stringent constraint is from the Belle experiment using 253 fb _1 of data. ^) 
Due to the neutrino, the final state is not fully reconstructed; however, at 
an e + e~ collider there are enough kinematic constraints to infer the neutrino 
momentum. Specifically, momentum conservation prescribes P v = P CM — 
P n Kl — P rost , where Pq M is the four-momentum of the e + e~ center-of-mass 
(CM) system, n s , K, and I are daughters from D* — > D°n s ^>Tr s K£i>, and P rest 
is the four-momentum of the remaining particles in the event. In the Belle 
analysis the magnitude |P rest | is rescaled to satisfy (P CM — P rcst ) 2 = , an d 
after this rescaling the direction of p rcst is adjusted to satisfy P 2 (= m 2 ) = 0. 

The AM = M v Klv - M Klv distributions for "right-sign" (RS) D° -» 
K~( + v and "wrong-sign" (WS) D — » K + l~v samples are shown in Fig. 1. 
Sensitivity to mixing is improved by utilizing information on the decay time, 
which is calculated by projecting the D° flight distance onto the (vertical) y 
axis: t = (M D0 /c) x (y vtx — Uip)/P y ■ This projection has superior decay time 
resolution, as the beam profile is only a few microns in y and thus the interaction 
point (y IP ) is well-determined. Events satisfying t>r D0 are divided into six t 
intervals, and the event yields and N$ s , acceptance ratio £yfs/ £ ns> anc ^ 
resulting mixing parameter are calculated separately for each. A^g and 
Nyj S are obtained from fitting the AM distributions. Doing a \ 2 fit to the six 
values gives an overall result r D = [0.20±0.47 (stat)±0.14 (syst)] x 1CT 3 , or 
r D — 0.10% at 90% C.L. No evidence for mixing is observed. The total number 
of signal candidates in all t intervals is 90601 ± 372 RS events and 10 ± 80 WS 
events. 
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Am [GeV] Am [GeV] 

Figure 1: AM distributions for RS D° — > K~i + u candidate decays (left) 
and WS D° — ► K + l~v candidate decays (right), from Belle using 253 fb~ x 

of data. ^) TTie W^S 1 pZoi shows no visible signal above background. 



3 7T+7T- CP-Eigenstate Decays 

When the final state is self-conjugate, e.g., K + K~, there is no strong phase 
difference between Af and Af. Assuming \Af\ — \Aj\ (no direct CPV), A = 
— \q/p\e 1 ^ and A = — \p/q\e~ l<t> , where is a weak phase difference and the 
leading minus sign is due to the phase convention CP\D°) = — \D°). Inserting 
these terms into Eqs. (1) and (2) and dropping the very small last term gives 



R(D°^K+K-) = \A 



K+K- 
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R(D°^K+K-) w |^ K+K _|2 e -rt e -b/ 9 |(ycos^+x S m«f^ ^ 

Eqs. (4) and (5) imply that the measured D° and D° inverse lifetimes are 
r[l + \q/p\(y cos 4>— xsm(f>)] and T[l + \p/q\(y cos <p + x sin </>)], respectively. We 
define y CP = T K -^+l T K + K - — 1, which equals \q/p\(y coscf) — xs'mcf)) for D° 
decays and \p/q\(y cos (f> + x sin (p) for D° decays. For \q/p\ = l, i.e., no CPV 
in mixing, y CP — ycoscj) for equal numbers of D° and D° decays together. If 
also = (no CPV), y CP = y. The observable y CP is measured by fitting the 
D° — > K + K~ and D° ^> K~ir + decay time distributions. 

To date, five experiments ^' ^' ^) have measured 2/ C p; the most precise 
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value is from BaBar using 91 fbr 1 of data. 8 ) To increase statistics, BaBar used 
both K + K~ and ir + ir~ decays, and, in addition, the D° — > K + K~ analysis used 
both a large inclusive D° sample and a smaller, higher purity sample in which 
the D° was required to originate from D* + — > D°tt + . The respective decay time 
distributions are shown in Fig. 2. Doing an unbinned maximum likelihood fit 
to each sample, combining results for K + K~ and 7r + 7r~, and taking the ratio 
of lifetimes gives y CP = [0.8 ± 0.4 (stat) Iq'I (syst)]%. This value is consistent 
with, but smaller than, the relatively large value measured by FOCUS: ^) 
y cp = [3.4 ± 1.4 (stat) ± 0.7 (syst)]%. 

BaBar also measures AY" = (r+ — t~)/(t+ + t~) x t k _ 7v+/ /{t), where 
r+ (r") is the lifetime for D°^K + R- (D°->K+K-) and (r) = (t + +t~)/2. 
For \q/p\ = l, AY = xsin<f>. The result is AY = [-0.8 ±0.6 (stat) ±0.2 (syst)]%, 
which indicates that either x is small or <p is small. 



4 D°{t) — »fC+7r Doubly-Cabibbo-Suppressed Decays 

For L>° K+tt-, A f is DCS, ^ is CF, and thus |^| < [.4^. In addi- 
tion, there may be a strong phase difference (S) between the amplitudes. 
Defining R D = \A f /A f \ 2 and R D = \A f /A f \ 2 , A = \q/p\ R D 1/2 e^ +<5 ) and 

X=\p/q\ R D e l (-<f>+ s ). Inserting these terms into Eqs. (1) and (2) gives 
R(D"^K+tt-) cx 
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,-rt 
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R D [y' cos (f) + x' sin <j>\ (Tt) 



(x 12 + y' 2 ) 



(nf 



V 2 + y' 2 )^2 



-(Ttf 



(6) 
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where x' = x cos 5 + y sin 5 and y' = —x sin S + y cos 8. These "rotated" mixing 
parameters absorb the unknown strong phase difference S. CPV enters Eqs. (6) 
and (7) in three ways: \q/p\ ^ 1 (CPV in mixing), R D ^ R D (CPV in the 
DCS amplitude), and <fr ^ (CPV via interference between the DCS and mixed 
amplitudes). Assuming no CPV gives the simpler expression 



R cx e 



-rt 



(x' 2 +y' 2 ) 



(Tt) 2 



(8) 



■5 




-4 -2 2 4 6 -4 -2 2 4 6 

Decay Time (ps) 



Figure 2: Decay time distributions for CF D — > A" ir + (upper left), D° — > 
K + K~ (upper right), D°^7r + 7r~ (lower left), andD°^K + K~ selected with- 
out using a D* + tag (lower right), from BaBar using 91 fb^ 1 of data. ^) TTie 
shaded histograms show the signal component obtained from the fit; residuals 
from the fit are plotted below each distribution. 

To date, six experiments ' ' 

12) 

have done a time-dependent analysis 
of D° — > K + ir~ decays; the most stringent constraints on x' 2 and y' are from 
Belle using 400 fb -1 of data. ^) The reconstructed M K7r and Q distributions 
after all selection criteria are shown in Fig. 3; fitting these distributions yields 
1073993 ± 1108 RS signal events and 4024 ± 88 WS signal events. Those 
events satisfying \M K7r — M D0 \ < 22 MeV/c 2 and \Q - 5.8 MeV| < 1.5 MeV 
(4cr intervals) have their decay times fitted for x' 2 , y', and R D . The results are 
listed in Table 1; projections of the fit are shown in Fig. 4(left). 

A 95% C.L. region in the x l2 -y' plane is obtained using a frequentist tech- 
nique based on "toy" Monte Carlo (MC) simulation. For points a — (x' 2 , y'), 
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Tabic 1: Limits on mixing parameters obtained from fitting the decay time 
distribution of WS D° — > K + tt~ decays, from Belle using 400 fiT 1 of data. ^) 



Fit Case 




Fit Result 

(xio- 3 ) 


95% C.L. interval 

(xio- 3 ) 




x' 2 


o.is±8:S 


< 0.72 


No CPV 


y' 


o.6lS:8 


(-9.9,6.8) 




R D 


3.64 ±0.17 


(3.3,4.0) 








(0.63 xl0~ 5 , 0.40) 




x 12 




< 0.72 




y' 




(-28,21) 


CPV allowed 






< 0.40 




A D 


23 ±47 


(-76,107) 




A M 


670 ± 1200 


(-995, 1000) 


No mixing/ CPV 


R D 


3.77 ± 0.08 (stat) ± 0.05 (syst) 



one generates ensembles of MC experiments and fits them using the same proce- 
dure as that used for the data. For each experiment, the difference in likelihood 
AL = lnL max — lnL(a) is calculated, where L max is evaluated for x' 2 > 0. The 
locus of points a for which 95% of the ensemble has AL less than that of the 
data is taken as the 95% C.L. contour. This contour is shown in Fig. 4(right); 
projections of the contour are listed in the right-most column of Table 1. 

CPV is accounted for by fitting the D° — > Kit and D — > Ktt samples 
separately; this yields six values: x' 2± , y^, and R^. Defining R M = (x' ±2 + 
y l±2 )/2 and A M = (R+ - R M )/(R M + R M ), one finds 



x 



I 1 ^ A M 
l±A 



1/4 

(x'cos^ ± y' sine))) (9) 



1/4 

y>± = {wa~) {y ' C0S(p T x ' s[n ^> ( 10 ) 

where there is an implicit sign ambiguity in x'^ due to Eqs. (6) and (7) being 
quadratic in x' . To allow for CPV , one obtains separate 1— V0. 05 = 77. 6% C.L. 
contours for (x' +2 , y' + ) and (x'~ 2 , y'~); points on the (x' +2 , y l+ ) contour are 
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then combined with points on the (x'~ 2 , y'~) contour and the combination 
used to solve Eqs. (9) and (10) for x' 2 and y'. Because the relative sign of x' + 
and x'~ is unknown, there are two solutions (one for each sign); Belle plots 
both in the (x 12 , y') plane and takes the outermost envelope of points as the 
95% C.L. contour allowing for CPV. This contour has a complicated shape [see 
Fig. 4(right)] due to the two solutions. Projections of the contour are listed in 
the right-most column of Table 1. In the case of no CPV, the no- mixing point 
x' 2 = y' = lies just outside the 95% C.L. contour; this point corresponds to 
3.9% C.L. with systematic uncertainty included. 
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Figure 3: WS D° — > K + tt decays: M K ^ spectrum for events satisfying 
Q G (5.3, 6.5) MeV (left), and Q spectrum for events satisfying M K7T S 
(1.845, 1.885) GeV/c 2 , from Belle using 400 fb' 1 of data. 12 ) 



5 D°(t)^K° tv+tv- Dalitz Plot Analysis 

In this method one considers a self-conjugate final state that is not a CP 
eigenstate, e.g., a three-body decay that can have either L — Q (CP-even) or 
L = l (CP-odd). If CPV is negligible, CP-eigenstates (denoted D_, D + ) are 
mass eigenstates (denoted D x , D 2 ), and the amplitude for D°(t)— >Kg 7r + 7r~ is: 

A K o^ = ±({K° s n+n-\H\D_(t)) + {K° s 7r+n-\H\D + (t))) 

= _4 e -(r 1 /2+im 1 ) t _|_ e -(T 2 /2+im 2 )t (--q-j 
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Figure 4: Left: WS D° — + K + ir decay-time distribution and fit projections. 
Right: 95% C.L. region for x' 2 , y' . From Belle using 400 fb^ 1 of data. ^) . 



=* R KOnv = \A_\ 2 e-^ 1 ^ + \A + \ 2 e- r ^ + 

2e~ Tt [Re(A + A*_)cos(Amt) + Im(.4 + .4* ) sin(Amt)] , (12) 

where A + _ is the amplitude for D + Kg tt + ir~ multiplied by l/(2p). Note 
that x= (m 2 — m 1 )/r and y= (T 2 — T 1 )/(2r). For a three-body final state, one 
can distinguish the A + and A_ components via a Dalitz plot analysis; i.e., a 
K s / (980) intermediate state is CP-even and contributes to A + , Kg p° is CP- 
odd and contributes to A_, K*(890) + ir~ is a flavor-eigenstate and contributes 
to both A + and A_, etc. Thus one models A + _ by separate sums of amplitudes 
5^ • aj e l5 i Aj , where Aj is the Breit-Wigner amplitude 

13) 

for resonance j and 

is a function of the Dalitz plot position M^.„ + , Mf COn _ . Using the probability 
density function of Eq. (12), one does an unbinned maximum likelihood fit to 
M|- 0jr+ , M^- 07r _, and the decay time t to determine a,j, Sj, x, and y. There 
is systematic uncertainty arising from the decay model, i.e., one must decide 
which intermediate states to include in the fit. Unlike Eq. (6), Eq. (12) depends 
linearly on x (x<l) and is therefore sensitive to its sign. 

This analysis was developed by CLEO, and their result 

14) 

based on 

9.0 fb _1 has not yet been superseded. To minimize backgrounds, the D° can- 



9 



didate is required to originate from D* + — ► D°7r + . The final Dalitz plot sample 
(Fig. 5) contains 5299 events with only (2.1 ± 1.5)% background. 

15) 

The decay model used consists of D Q -» K*(890)~ n+ K*(U30)o 2 tt+ 
Jf*(1680)-7r+ ^/ (980), tf°/ 2 (1270), tf°/ (1370), WS D° - 

if *(890) + 7T _ , and a nonresonant component. The fit results are listed in Ta- 
ble 2; the 95% C.L. intervals correspond to the values at which — 21n£ rises by 
3.84 units, where C is the likelihood function. CPV is included in the fit by 
introducing parameters e = (p— q)/(p + q) (in analogy with K° decays) and </>, 
the weak phase difference between A K o V7r and A K o nn - The results listed are 
consistent with no mixing or CPV. 
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Figure 5: Dalitz plot (lower right) and projections (lower left, upper plots) for 
D°^K%ir+n- decays, from CLEO using 9.0 ft- 1 of data. 15 ) 
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Table 2: Limits on mixing and CPV parameters from a t-dependent fit to the 
D° -> K^tt+tt- Dalitz plot, from CLEO using 9.0 fir 1 . 1J ^> The errors are 
statistical, experimental systematic, and modeling systematic, respectively. 



Fit 


Param. Fit Result (%) 95% C.L. Inter. (%) 


No CPV 


x 1.8 ±H ±0.4 ±0.4 (-4.7,8.6) 
y -1.4^ ±0.8 ±0.4 (-6.3,3.7) 


CPV 
Allowed 


x 2.3 ttl ±0.4 ±0.4 (-4.5,9.3) 
y -1.5 tH ±0.8 ±0.4 (-6.4,3.6) 
e 1.1 ± 0.7 ± 0.4 ± 0.2 (-0.4, 2.4) 
4> (5.7 ±2.8 ±0.4 ± 1.2)° (-0.3°, 11.7°) 



6 D°(t)^K+Tv-n° and K+n-n+n- Multibody Decays 

Mixing has also been searched for in WS multibody final states u ' 

16, 17) 

K + TT~ir a and K + ir~Tr + ir~; the most recent measurement is from Belle using 
281 fb" 1 of data. 17 ) The final signal yields arc 1978 ± 104 D° -> K+ir-ir 
decays and 1721 ± 75 D° — > K + Tr~-K + iT~ decays. For this analysis no decay 
time information is used, i.e., Belle measures the time- integrated ratio of WS 
to RS decays: 

Kws ~ jR[D^K-ir+(n^)]dt ~ H d + V K dV + . W 

where R D is the ratio of the DCS rate to the CF rate as previously defined for 
D°^K + n- decays. The results are R ws = [0.229 ±0.015 (stat) ±S:8oi (syst)]% 
for K+ir-ir and [0.320 ±0.018 (stat) ±g gi| (syst)]% for K+tt-tt+tt-. Inserting 
these values into Eq. (13) allows one to determine R D as a function of x' or y'. 
Assuming x' = and \x'\ =0.027 gives the curves shown in Fig. 6; the latter |x'| 
value corresponds to Belle's 95% C.L. upper limit from D° — > K + tt~ decays 
(see Table 1). However, the value of x' from D° — > K + ir~ may differ from that 
from D° — > K + ir~n(n) due to the strong phase differences (S) being different. 

7 Summary 

The 95% C.L. allowed ranges for x' and y' are plotted in Fig. 7; for simplicity 
we assume negligible CPV. The most stringent constraints are \x'\ < 2.7% 



11 



0.8 




-0.06 -0.04 -0.02 0.02 0.04 0.06 



y 

Figure 6: WS D° — > K + n~n(ir) decays: 95% C.L. range for R D as a function 

ofy' for \x'\ = and \x'\ = 0.027, from Belle using 281 ft)- 1 of data. V/ ) The 
point with la error bars is Belle's result from D° — > K + n~ decays (see Table 1). 

and y' e (—1.0%, 0.7%). These ranges are projections of the two-dimensional 
95% C.L. region for x' 2 , y' from Belle [Fig. 4(right)]. 

The results for y^p are plotted in Fig. 8. Here the central values and 
lcr errors are shown; combining the results assuming the errors uncorrelated 
gives y CP = (1.09 ± 0.46)%. This value differs from zero by 2.4cr and indicates 
a nonzero decay width difference Ar. Assuming negligible CPV, one can 
combine this value with Belle's central value for y' , (0.06 io' 3 g)%- The result is 
y'/y = cos 5 — (x/y) sin 5 — 0.05 ^37, wnere the- error is obtained from an MC 
calculation as the fractional errors on y and y' are large. This small central 
value (albeit with a large error) implies tan S ss y/x; i.e., if x <C y, then 5 is 
near 90°. Such a strong phase difference would be much larger than expected. 
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